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ABSTRACT 

We have examined In detail the synchrotron emission by electrons of 
energy greater than a few TeV in Earth's magnetic field. The photon spectrum 
lies in the X-ray and y-ray region. As the emission takes place in a narrow 
cone along the direction of electron, the photons would be incident nearly 
along a straight line on a detector. This unique feature provides the 
signature to identify the electron unambiguously. The mean energy of the 
photons being proportional to the square of the electron energy allows us to 
determine the energy accurately. Though it may appear that one needs to know 
the arrival direction of electrons to obtain its energy, we have shown that an 
omnidirectional detector can be satisfactorily used to estimate the energy. 

We also show that the collecting power of the detector is a sensitive function 
of the area of the detector A, the energy of electron Eq, and the number of 
photons required to identify an electron n^; asymptotically the collecting 
power Is proportional to A E^n” . An instrument, with a energy 
threshold for the detection of photons, can be used to measure reliably the 
integral flux of electrons, even if it has limited energy resolution. We have 
calculated the event rate expected using an ideal balloon-borne detector 
capable of detecting above 20 keV at 4 g cm”^ of atmospheric depth over 
Palestine, Texas, and compared with the expected rates using instruments based 
on currently available technique of detection. 
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1. INTRODUCTION 

High energy cosmic ray electrons play an Important role In the 
understanding of the propagation of cosmic rays. This stems from the fact 
that, unlike nuclei, electrons have a very low mass and as a result, they 
suffer severe energy loss especially at high energies through synchrotron 
radiation and Inverse Compton scattering processes. The lifetime of 1 TeV 
electron against such energy loss processes. In a typical Interstellar 
environment. Is about a few times 10^ years; this lifetime decreases linearly 
with energy. As a result, the flux of high energy electrons are depleted. 

The shape of the energy spectrum of electrons at energies beyond a few TeV 
would be dependent upon the distribution of sources In the solar neighborhood 
and perhaps on the local dynamics of the propagation of cosmic rays. The 
Influence of a strong nearby source would also be seen as a marked enhancement 
over the depleted electron spectrum. This Information Is difficult to obtain 
from other components of cosmic rays and thus the detection of ultra high 
energy cosmic ray electrons has profound implications on our understanding of 
the propagation of cosmic rays In the solar neighborhood. 

A variety of experimental techniques has been employed to measure the 
energy spectrum of cosmic ray electrons below a few hundred GeV (see 
Balasubrahmanyan and McDonald, 1974). However, beyond this energy only large 
area nuclear emulsion detectors (Anand, Daniel, and Stephens, 1973; Nishimura 
et al., 1980) and a large area calorimeter (Muller and Meyer, 1973) have been 
made use of to extend the energy up to about 1 TeV. The limited collecting 
power of these detectors does not permit the extension of the electron 
spectrum beyond 1 TeV. A transition radiation detector. In combination with 
an Ionization calorimeter, placed on the space shuttle (Ormes, 1982) could 
reliably measure the energy spectrum only to about a few TeV. The limitation. 
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of this technique, to extend the measurements to higher energies comes from 
the fact that one needs a deep calorimeter and a transition radiator with 
sufficiently large spacings between elements (Cherry, 1978). The net results 
of these modification are to increase the total weight of detectors for the 
same area and to decrease the available solid angle. Thus, it appears that 
direct measurement of the energy spectrum of high energy electrons beyond a 
few TeV seems difficult with the existing techniques. 

We have also examined the possibility of detecting cosmic ray electrons 
by their shower development in the atmosphere; here one cannot distinguish 
between gamma rays and electrons. For a typical mountain altitude, one can 
make use of this technique for energies beyond 10^*^ eV. In order to eliminate 
air showers simulated by nuclear active particles to one part in 10®, one 
requires very large area detectors to monitor hadrons and muons associated 
with a shower. The required area of these detectors is such that the steady 
hadron and muon flux, resulting from the interaction of cosmic rays at all 
energies, is sufficiently large enough to be in coincidence with an electron 
event. As a result, the rare electron air shower would be mistaken for a 
nucleon initiated cascade and get vetoed. It is quite clear that in spite of 
the difficulty in making very large area detector arrays, it may not be 
possible to detect electron initiated showers in the atmosphere. 

High energy electrons lose energy through synchrotron radiation in the 
Earth's magnetic field. An electron of energy 10^^ eV radiates photons in the 
energy region between a few keV and a few MeV. These photons are emitted in 
the direction of the electron within an extremely narrow cone; they are 
distributed in the plane of curvature of the electron on the convex side of 
the curvature. As a result, these photons would be incident on a detector 
along a straight line with very little spread. This pattern provides a unique 
signature for the identification of the electron. Because of the fact that 


Identification of such a signature for an electron event does not require the 
electron to be incident on the detector, the response area of the detector 
could be larger than the detector area itself. This property provides large 
collecting power for the detector, which is very necessary when we look for 
high energy electrons. Further, the mean energy of the emitted photons is 
proportional to the square of the energy' of the electron. As a consequence, 
the energy estimation becomes reliable even if a detector has low energy 
resolution for photons. Thus, the detection of high energy electrons using 
this technique would open up an important window of energy for the study of 
primary cosmic ray electrons. 

In this paper we study in some detail the synchrotron emission by a high 
energy electron in the geomagnetic field and its dependence upon different 
arrival directions over Palestine, Texas, where major balloon-borne 
experiments are being conducted. We discuss how the response of a detector 
depends upon the arrival direction of electrons, the different criteria one 
adopts to identify an electron event, the area of the detector, and the energy 
of the electron. We then make use of an omnidirectional circular detector to 
examine whether it is possible to determine the energy of an electron without 
knowing its arrival direction. We are able to demonstrate the reliability of 
energy estimation using such a detector system. We then make estimates of the 
collecting power of a detector as a function of the energy of electrons for 
different detector areas with different selection criteria. Making use of 
this information we calculate the event rates expected by folding in the 
energy spectrum of cosmic ray electrons to show the viability of detecting 
cosmic ray electrons at energies greater than a few TeV. It may be pointed 
out that the performance of such a detector would have an advantage In the 
Equatorial region due to a favorable field configuration. 
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SYNC HROTRON RADIATION OF ELECTRONS IN EARTH^S MAGNETIC FIELD 
The theory of synchrotron radiation is extensively discussed in the 
literature, and here we follow the work of Ginzburg and Syrovatskil (1965). 

We assume in our treatment that the synchrotron results based on classical 
electrodynamics Is valid In the energy region > 10^3 ev with magnetic field 
strength of the order of a few tenths of a Gauss. The power radiated by an 
electron of energy EgiTeV) In a magnetic field strength B (Gauss) is 

P(Eq) = 3.77 X 10^ E^ MeV s‘^ (1) 

where B^^ 1s the component of the magnetic field perpendicular to the direction 
of electron. This radiation is peaked around an energy 

== 1.9 X 10"^ Bj^ Eq MeV. (2) 

If one assumes that most of the energy Is radiated around this energy, It Is 
clear from equations (1) and (2) that the mean number of photons emitted is 
proportional to the field strength B^ and is Independent of the energy of the 
electron. It Is, therefore, essential to look for this radiation in a region 
of high magnetic field strength. In the case of Earth's magnetic field with a 
typical field strength of Bj^ = 0.3 Gauss, the number of photons emitted per 
15 km of the electron trajectory is about 3 photons which is just sufficient 
to detect the presence of an ultra high energy electron. The energy of the 
photons emitted by a few TeV electrons Is in the X- and -y-ray regime, which 
can be detected by instruments used for the study of X- and Y-ray astronomy. 
The spectral shape of photons emitted by an electron can be written as 
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(E/E^) 


^ 5 / 3 ^x)dx photon/(MeV.s) (3) 


where “ 6.56 x 10“^ E^ MeV and x i; “ ^5/3< x)dx gives the spectral 
distribution of the power emitted by the electron; Is the Bessel 


function of the second kind. 


We have shown 1n Figure 1 the differential and Integral spectra of 
photons emitted by an electron In a magnetic field of strength B^ = 0.4G for 
electron energies 5, 10, and 50 TeV. A common feature that can be noticed 
from *he shapes of the differential spectra Is that at low energies, they have 
a constant slope of -2/3. The spectrum slowly steepens to a slope of -1,0 
around eIJ and rapidly drops off exponentially beyond Further one notices 
that the absolute flux of low energy photons decreases with Increasing energy 
of the electrons (curves B and C). This property Is useful In selecting a 
suitable detector In order that only a small fraction of the total radiated 
photons is undetected due to detector threshold. This is well demonstrated by 
the Integral spectra (dashed curves) shown 1n Figure 1, One notices from the 
dashed curve A that at 5 TeV more than 75 percent of the radiation takes place 
above 10 keV and only 5 percent above 1 MeV. This Implies that 70 percent of 
the photons emitted are spread over 2 decades 1n energy and 50 percent over a 
factor of 25 In the energy scale. It 1s quite evident from the dashed curves 
that they are nearly Identical, and they appear to be shifted In the energy 
scale by an amount proportional to the square of the electron energy. This 
points to the fact that. In order to cover the m decades of electron energy, 
one needs an Instrument capable of detecting photons over (m^ + 2) decades In 


energy. 

Synchrotron photons are emitted within a narrow cone In the direction of 
electron. The radiation is polarized, and the angular distributions of 
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photons along the two principal directions of polarization are different. The 
angular dependence of photons In the two principal directions of polarization 
Is written as (Ginzburg and Syrovatskll, 1965) 



(9) . (1 + [K 2/3 

(4a) 

and 

(e) « (e/*)^ (1 + (G^)]^ 

(4b) 

where 

^ (1 . e2/t2)3/2. 

(5) 

2 

In these equations, ij; = me /E and K's are the Bessel 

functions of the second 


kind with Imaginary arguments. The first component represented by equation 
(4a) Is polarized In the direction perpendicular to the projection of the 
magnetic field on a plane normal to the direction of motion of the electron. 
The second component represented by equation (4b) Is along the projection of 
the magnetic field. It Is clear from these equations that the first component 
peaks at 0 - 0 and drops off at larger values of e, while the second component 
does not contribute at 9 = 0 but peaks at a finite value of 0 . The relative 
contribution of these two components Is dependent upon the photon energy and 
the magnetic field strength through equation (5) and through the dependence of 
Ej, on Bj^. One may point out that the angle of emission In Inversely 
proportional to the energy of the electron. 

In our calculations, we first determine the number of photons In finite 
energy intervals through equation (3) and then calculate their angular 
distributions. Therefore, one needs to obtain the normalized angular 
distributions of the photons for various values of This Is obtained by 

Integrating equations (4a) and (4b) as 
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The normalized angular distribution of the radiated photons are shown In 

Figure 2 as a function of (9/i{;) for various values of (E^/Ec^’ clear 

from this figure that the angular distribution is a sensitive function of the 

energy of the photon; the smaller the value of E^/E^* the wider is the 

distribution. This behavior results from the fact that the second component 

becomes important (equation 4b) for small values of E,^, as it is evident from 

the appearance of a broad peak for 't /E = 0.05. Because of the fact that 

most of the energy radiated is peaked around E /E “ 0.29, the wider 

Y c 

distributions, which are characteristic of the low energy photon emission, do 
not contribute to the integral photon distribution. This is demonstrated in 
Figure 3, in which we have plotted the angular distribution of photons of 
energy > 20 keV emitted by an electron of energy Eq = 10 TeV in a magnetic 
field = 0.24 Gauss. Curve A in Figure 3 is the number of photons of energy 

> 20 keV emitted per unit radian by an electron over a path length of 10 km as 
a function of angle from the direction of motion of the electron. In this case, 
the value of E^. = 1.57 MeV and the threshold energy of 20 keV corresponds 
to E^/E^ = 0,013; the fraction of photons emitted below 20 keV is about 15 
percent. One notices that at small values of 0, the angular distribution of 
all photons resemble the distributions corresponding to large values of 
E /E in Figure 2 and becomes wider as e increases, due to the contribution 

Y C 

from lower energies. The integral photon distribution N^(< e) within a cone 
angle 0 is shown by Curve C in Figure 3. From this curve one infers that 70 
percent of the photons of E^ > 20 keV are emitted within a narrow cone of 0.65 
X 10"^ radian. This correspond to a spread of only 1.3 cms on to a detector 
placed perpendicular to the path of an electron, if the radiation is emitted 
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from a distance of 100 km. Similarly, 90 percent of the radiation falls 
within 2.3 cms. 

The angular distribution is narrower for higher energy photons. As a 
result, most of the radiated energy is concentrated in a small cone. We have 
shown by Curve B the energy flow I E^(0) per unit radian as a function of e. 
From this curve, one notices that 90 percent of the energy is contained in a 
narrow cone of 4.5 x 10’^ radian, which corresponds to a spread of 0.9 cm onto 
a detector kept at a distance of 100 km from the point of emission. The 
emission cone decreases linearly as the energy of electrons increases. The 
effect of increasing the magnetic field increases both the number of photons 
and their energies. Curve A" In Figure 3 shows the angular distribution when 
the magnetic field is increased by a factor of two (0.48 Gauss). Since E^/Eq 
is unaffected by the increase of the magnetic field, the angular distribution 
is nearly unaffected. The flattening of the distribution at large e values 
result from the inclusion of photons of energy very much lower than E^, because 
these photons are now above the low energy threshold set at 20 keV. 

From the above discussions it is clear that the spread in the photon 

distribution is small, compared to the required size of the detector, and one 

can consider the photons to be distributed approximately along a straight line 
in the detector. This also Implies that very high spatial resolution is not 
needed to detect electrons. However, higher resolution would provide 
information on the photon density and on the linear spread, both of which can 
be made use of in the determination of the energy of the electron. The mean 
energy of the photons emitted by the electron can be calculated using Equation 
(3), and it is found that 

<E > = 0.02 B , E^ MeV 

Y 10 


( 7 ) 
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which Is approximately the energy at which the power radiated by the electron 
is peaked (equation Z), If the detector has a low energy threshold, the mean 
energy would be larger; high energy cut-off would give a lower value, 

3* ‘MOTION OF ELECTRONS IN EARTH^S MAGNETIC FIELD 

In this section we discuss the method of tracing the trajectory of cosmic 
ray electrons In the Earth's magnetic field, We describe briefly the 
procedure adopted to calculate the synchrotron emission by the electron as it 
traverses the geomagnetic field, Ths equation of motion of a charged particle 
in a magnetic field is given by 

W = Sc <3> 

Where e, m and v are respectively the charge, mass, and velocity of the parti cl 
This equation can be converted to a system of six simultaneous linear 
differential equations in spherical coordinates (McCraken et al, 1962) so that 
one can numerically integrate these equations to trace the trajectory. For 
the magnetic field term 5 in Equation (8), we have made use of the sixth 
degree simulation of the geomagnetic field with the field coefficients given 
by Finch and Leaton (1957). This is done in order to take advantage of the 
computer program of McCraken et al. (1962) to determine the components (v x 5) 
in the spherical coordinate system using the complex geomagnetic field 
potential (we thank Dr. M. A. Shea for this information). It may be mentioned 
here that, at energies of interest here (> 10^^ eV), the use of different 
geomagnetic field coefficients would not alter any of the basic conclusions 
derived in this paper. 

We have chosen a typical balloon altitude corresponding to a depth of 
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^ g of atmosphere over Palestine, Texas as the initial point on the 

trajectory of the electron. The detector plane is assumed to be horizontal. 
Trajectories with different arrival directions (zenith Z and azimuth East of 
North are traced from the detector plane outwards. For this purpose, 
stepwise integration was carried out similar to the work of Daniel and 
Stephens (1966). We have shown in Figure 4a an exaggerated view of the 
trajectory of an electron. In this figure 0 is the arrival point, B's denote 
the steps with step length As, and A's are the intercepts on the detector 
plane by tangents BgAg,... at point B's on the trajectory. The wavy 
lines are indicative of the photons emitted by the electron. These photons 
intercept the detector within s. narrow angle 6. At each step, the following 
aspects are examined. 

a. Ffdf the new coordinates, the local vertical height and the 
angular deviation are obtained. 

b. The total magnetic field strength B and its component 
perpendicular to the direction of electron B^^ are determined. In our 
location, the field is directed approximately towards 30® zenith and about 
190® azimuth, Using the value of B^^, the radius of curvature Rj,(S) at S 

- I S-AS/2, the distance swept by the tangent ah = tan {S/Rj,(S)} 
x AS/cosZ, and a = AH (see Figure 4a) are calculated. The angular deviation 
obtained using procedure (a) is checked against tan“^ {As/R (s)}. 

c. The energy spectrum of synchrotron photons is obtained using Eq 
and Bj^. The attenuation length of 20 keV photon is about 3 g cm"^ of air and 
this decreases rapidly with decreasing energy. Therefore, we have considered 
only photons of energy greater than 20 keV in this analysis. The photon 
spectrum above 20 keV is then corrected for the attenuation in the air between 
the point of emission and the detector. The matter traversed by the photons is 
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obtained from the standard atmosphere using the relation x = exp{(48.5-h)/7.54} , 
where x Is the atmosphere depth In g cm"^ and h Is the vertical height In km; 
this relationship Is valid for x < 11 g cm"^. We have made use of the local 
zenith angle to correct for the slant depth and the absorption cross-sections 
compiled by Zombeck (1980) to correct for the absorption. 

d. The total number of photons received over AA, the total energy 
flow, the mean energy <E^>, and their angular distribution are calculated. 

Using these, we obtained the quantities N^{< a) and E^(< i). These are found 
to be proportional to where 3 = 0.5 for small values of A, decreasing 
slowly with increasing A due to the decrease In the magnetic field strength. 

e. Finally, the linear density of photons over aa Is calculated, and 
the Integration Is continued until this density falls below a preset value. 

With the above Information we have proceeded to calculate the response of 
a detector to determine the specific characteristics of a detector. 

4. RESPONSE OF A DETECTO R 

It can be noticed from Figure 4a that, when an electron passes through 
the detector, all the photons would be Incident only on the convex side of the 
trajectory along a straight line. This Is Indeed a special situation and 
makes the Identification of the electron event very unique. If the electron 
does not intercept the detector, one can still make use of the property that 
the Incident photons are col 11 near In the detector to detect electrons. The 
number of photons required to Identify an event unambiguously Is n^ = 3. 

However, depending upon the resolution of the detector, one may need a 
different criterion with n^ > 3. This criterion sets constraints on the size 
and spatial resolution of the detector. Therefore, we examine In some detail 
the response of a detector as a function of energy for different values of n^ 
and for various detector sizes. 
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In order to Investigate the response of a detector, we define the 
response area of an Instrument as the area in the plane of the detector, in 
which the instrument detects at least n^ photons emitted by the electron 
arriving within this area. By virtue of this definition, it is clear that Af^ 
Is not the same as the detector area A, but could be smaller or larger. The 
following properties of Af^ can be inferred from the theory of synchrotron 
radiation and from Figure 4a, in which one notices that the linear density of 
photon on the detector plane decreases as the distance A from the electron 
increases. 

a. As the value of n^ increases, one requires the detector to be 
close to the incident electron and hence, Ar decreases with Increasing 

b. The bigger the detector, the higher is the number of photons 

detected even if the detector is kept at a large distance from the electron, 

and therefore increases with A, 

c. The mean number of photons emitted by an electron is proportional 

to B , while the radius of curvature is inversely proportional to B . As a 

result, the linear density of photons emitted (Figure 4a) is independent of 
Bj^. As the displacement a of the electron increases with B^^, Ar increases 
with B^ hence depends on the arrival direction of electron (Z ,(}>). 

d. As the energy of the electron increases, the curvature decreases 
and the linear density of photon increases. This enables higher energy 
electrons to be detected even at large a from the detector, and hence Ap 
depends upon Eq. 

e. The photon spectrum decreases exponentially above (Figure 1), 
Therefore, for an instrument with a threshold energy for detecting photons, Ap 
would decrease to zero rapidly below a critical low value of Eq, at which the 
threshold energy becomes larger than Eq. 
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From the above, It is clear that is a function of A, n^, 1,<P, and 
Fq. In order to evaluate these dependences quantitatively, we have made use 
of a circular detector so that the circular symmetry eliminates the influence 
of shape. As shown in Figure 4b, the detector is divided into a number of 
strips of width Ab parallel to the plane of curvature of the electron. Let 
the length of one such strip be L and n^ = 3. If the electron arrives at 0, 
it would have 3 photons (filled circles) associated with it to be detected by 
the instrument between 0 and D. Since no photons would be incident on the 
concave side of the trajectory, part of the detector stip OA does not record 
any photon. Therefore, the minimum length required by the detector is OD = 

The value of Ji, can be determined from N^(< distribution, described in 
section 3d, by setting n^ = (< z^). As the electron moves towards A, the 

number of photons detected increases, reaching a value maximum at A, and then 
decreases. At O' the number of photons recorded by the entire strip is 3 
(open circles). The electron arriving beyond O' cannot give rise to 3 photons 
to be detected and, hence, the maximum distance up to which an electron 
can be detected is O'D = z^. The value of Z 2 is determined from the 

relation n = N (< z^) - N {< iz^- L)}. The response area of the strip is 

y y u y c 

(Z 2 ~ X Ab, which is different from the area of the strip. The total 
response area is the sum of the response area determined for all paralled 
strips in the direction of (p. 

We have illustrated in Figure 5 the response pattern of a circular 
detector of radius 150 cms for electrons arriving at a zenith angle of 45° 

over Palestine, Texas. Figure 4a is for Eq = 10 TeV. The detector is shown 

by thick circles. The thin curves trace the circumference of the response 
area for n^ =5 in the direction of the North, East, and West. The response 

area towards the South is zero because the magnetic field direction is only 
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15® from the arrival direction and hence is small, The effect of 
increasing the value of n^ to 7 photons is shown by the dashed curve for 
North. The dramatic increase of the response area with Eq is shown in Figure 
5b for Eg = 20 TeV; here for convenience, this figure is scaled down by a 
factor of 2, As in Figure 5a, the dotted curve represents the response area 
of the detector for = 7. Note that at this energy there is a finite 
response area in the direction of South, because as Eg increases, the energy 
of the photons increases above the threshold set at 20 keV. The following 
inference can be obtained from this figure; The response area (a) varies with 
the arrival direction of electrons, (b) appears elongated along the diameter, 
indicating its dependence on the shape and size of the detector, (c) decreases 
as the value of n^ is increased, and (d) increases with energy. 

The dependence of Aj^ on the arrival direction of electrons is shown in 
the upper part of Figure 6 for the same circular detector and for Eg = 

10 TeV. Curves A, B, C, and D correspond to zenith angles 30®, 45®, 60®, and 
75®, respectively for n^ = 5 photons; curve A' is for n^ = 7 photons. In the 
middle part of the figure, we have plotted value as a function of the 
azimuthal angle. One notices a general similarity between these curves and 
those of Aj^. However, the absolute values of A|^ at large zenith angles are 
smaller than at small zenith angles, unlike the values of Bj^. This is due 
to the fact that the projected area of the detector perpendicular to the 
direction of electron decreases with cosZ. If one takes this into account, A|^ 
curves would resemble a plot of (B^- e) as a function of t|), where e is a 
constant. Thus one can represent the response area as, 

Ap(Z, = a[B^ (Z, ^iCosZ - bCos'^Z]^ 


(9) 
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Here the constants a and b are function of Eq, n^, and A; a Is approximately a 
reciprocal of n^ and Is weakly dependent upon Eq, and 3 Is nearly an 
Independent constant. For the example given In Figure 5, the values of the 
constants are a = 27 m^, b = 0.122, a = 0.2, and 3 = 0.75, 

In the lower part of Figure 6 we have plotted the mean number of photons 
> 20 keV emitted by a 10 TeV electron over a path of 10km. These curves 
resemble the variation of Note that the photon emission rate is rather 
small. Therefore, one needs a sufficiently large area detector for an 
experimental Investigation of high energy electrons, apart from the fact that 
the flux of these particles Is very small. 

The response area of a detector as a function of zenith angle is plotted 
In Figure 7, for various azimuthal angles and for different primary 
energies. The azimuthal angles chosen In this figure are North, East, 

Southeast, South, Southwest, and West, which are shown by curves A, B. C, D. 

E, and F, respectively. The top part of this figure Is for primary energy Eq 
= 10 TeV, the middle part for 5 TeV and the bottom part for 30 TeV; all these 
curves are for = 5 photons for a detector with radius 150 cms. It Is clear 
from this figure that as the primary energy Increases, the response area extends 
to arrival directions at large zenith angles and towards South (low B^). It 
Is suggestive from this figure that If one sets zero level for A[^ at around 

9 m^ for 30 TeV, the shapes of these modified curves would resemble those of 

10 TeV. Similarly a zero-level setting at 6m^ for 10 TeV would make these 
curves resemble the 5 TeV curves. These zero-level settings are the result of 
detector threshold, by which low energy photons emitted by the electron, 
either due to Eq being small or due to low B^^, are not detected. Below 5 TeV, 
the response area decreases rapidly and becomes negligible even at 4 TeV. 
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This is an important characteristic of a detector with low energy threshold 
for photon detection. 

We have demonstrated in this section how the response of an instrument 
depends upon the arrival direction of the electron, the minimum number of 
photons required to identify an electron event, the active area of the 
detector, and the energy of the electron. It is clear from our discussions 
that for a given Eq there is a threshold area for an instrument, below which 
it cannot detect electrons. We have also shown that a detector with a low 
energy threshold for detecting photons, provides a natural threshold energy 
for electrons. One can take advantage of this property to measure integral 
flux of electrons. We may also point out that the response area of an 
instrument depends also on the shape of the detector, as is evident from 
Figure 5. It can be shown that the response area is larger for a detector 
elongated in the East-West direction than for a circular detector with the 
same area. However, the advantage of an elongated detector needs careful 
study because of the difficulty in design and orientation. 

5. ENERGY ESTIMATION OF ELECTRONS 

The mean energy of photons emitted by an electron (Equation 7) is 

2 

proportional to and to E^. Therefore, for a given the energy of an 
electron can be determined from the observed mean energy of photons. It is 
clear from Figure 6 that the value of varies with the arrival direction, 
indicating the necessity of determining the direction of photons. On the 
other hand, in the X- and y-ray regime under consideration, it is difficult to 
build large area detectors, capable of providing directional information with 
good detection efficiency. However, one also infers from the same figure that 
the response area increases with the increasing value of B^^. The latter 
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observation suggests that the major fraction of the electrons detected by the 
instrument is dominated by arrival directions with large B values. This 

4 * 

interesting property helps us to examine the viability of using an omnidirec- 
tional detector to study cosmic ray electrons. 

The number of particles detected by an instrument in a direction (Z, •}>) 
is given by 

MEq, Z, = CP{E^, Z. <^)tj(E^)dE^ (10) 

where j{£^) is the differential energy spectrum of particles, t is the time of 
observation and CP is the collecting power of the instrument. The collecting 
power is the product of the projected response area perpendicular to the 
arrival direction and the solid angle is that direction. Since the number of 
particles observed is directly proportional to CP, it is possible to draw 
histograms of the calculated mean energy <E^> using CP. For this purpose, we 
have divided the zenith angles in to 5® intervals and the azimuth angles into 
15® intervals, and calculated for each solid angle bin the values of CP and 
<E^>. Using these values we have constructed histograms of <E^>. 

We have plotted in Figure 8 histograms of <E^> using an omnidirectional 
detector of radius 150 cms for Eq = 10 TeV and n^ = 5 photons. The 
distributions marked as A, B, C, and D correspond to arrival directions with 
zenith angles less than 30®, 45®, 60® and 85®, respectively; they include all 
azimuthal angles. Because of the finite size of the angular intervals chosen 
for this purpose these distributions do not look smooth. One notices from 
this figure that the peak in the distribution shifts to higher energies as one 
includes larc;er zenith angles. This is because the mean value of B^^ over all 
azimuthal angles is higher for large zenith angles than for small Z (see 
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Figure 6). A closer examination of Figure 8 reveals that the dispersion about 
the peak Is nearly the same for all distributions. They are asymmetric with 
half widths at half maximum of about 27 percent of the peak value towards 
lower energies and about 10 percent towards higher energies. The sharp cutoff 
In the distribution at higher energies is an important feature because while 
determining the energy spectrum, It minimizes the spillover from the large 
flux of low energy electrons to higher energies. 

The distributions of mean energy In Figure 8 shows that the energy 
corresponding to the peak In the distribution depends upon the acceptance 
angle of the Instrument. It Is also Important to Inquire w'^^ther this peak 
energy depends upon other parameters such as n^ and A. For this purpose we 
have made use of the calculated collecting powers for zenith angles less than 
60®. <E^> distributions for various values of Oy are shown in Figure 9; the 

other parameters are kept the same as In Figure 8. The distributions marked 
as A, B, and C correspond to n^ values 3, 5, and 7 photons, respectively. It 
Is clear from this figure that the peak value of the energy is the same for 
all the distributions, suggesting that the peak energy is Independent of n^. 
Note that the width of the distribution Is also Independent of ny A closer 
examination of this figure reveals that, while the distribution terminates 
sharply at higher energies for all values of n^. It extends to lower energies 
as r\y decreases. The former observation Is the result of the zenith angle 
restriction, and the latter 1s due to the inclusion of low directions. 

We have shown In Figure 10, the histograms of <Ey> for various detector 
sizes. Histograms marked as A, B, anc C correspond to circular detectors of 
radii 100 cms, 150 cms, and 200 cms, respectively. We have restricted the 
zenith angles to < 60® and kept n^ = 5 photons, so that the histograms C in 
Figure 8, B In Figure 9, and B in Figure 10 are the same. Figure 10 indicates 
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that the peak energy as well as the width of the distribution are independent 
of the area. We also notice that as we increase the area beyond a certain 
limit, the <E^> does not extend to lower energies, Thus, a comparison between 
Figures 9 and 10 suggests that the extent of the distribution to lower 
energies is a sensitive function of n^ rather than of A. 

It is quite evident from the above description that one can determine for 
a detector, the most probable value of <E^> using the calculated histogram of 
the type described above. Accordingly, we have estimated for a balloon-borne 
instrument with 20 keV cutoff, the most probable value of the mean energy 
<£y>m at a float altitude of A g cm"^ of atmospheric depth over Palestine, 

Texas. A plot of the estimated <EY>m as a function of the electron energy is 
shown in Figure 11, It is seen that the calculated data points can be 
represented by a power law of the type, 

<E^>j^ = 0.0225 eJ’^® MeV (11) 

As inferred earlier, this relationship should be independent of n^ and A. It 

can be seen from Figure 11 that at low values of Eq below 5 TeV, the 

calculated value deviates from the power law as a result of the low energy 

threshold of the detector and atmospheric absorption of X-rays. The energy 

dependence obtained here is flatter than the one given in Equation (7). This 

is because as the energy increases, the detector responds to photons emitted m 

fi 

from greater distances from the detector, where the magnetic field is ji 

}• 

weaker. It is Important to mention here that the above relation in Equation 
(11) is based on the assumption that the instrument has 100 percent efficiency |i 

to detect and to measure the energy of photons > 20 keV, without high energy | 

cutoff. In actual practice this is not true, and one needs to obtain a | 

I 
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similar relation by including both the detector efficiency and the energy 
response. 

In an experiment, one makes use of a relation similar to equation (11) to 
convert the observed mean energy of the photons to electron energy. In the 
same manner, we have converted the calculated mean energy in each angular bin 
to electron energy and examined the resultant distribution using the same 
procedure described earlier. We have plotted in Figure 12, the collecting 
power as a function of the ratio of the estimated energy to the actual 
•electron energy (E/Eq) in the form of a histogram, for different Eq values. 
Histograms marked as A, B, C, and D correspond to electron energies 5 TeV, 

10 TeV, 20 TeV, and 50 TeV, respectively. The value of n^ chosen for this 
plot is 5 photons and the radius of the detector is 150 cms; no angular 
restriction is introduced in this plot. 

From Figure 12, one infers that except at low primary energies, where the 
threshold effect narrows the distribution, the rest of the distributions are 
similar. The apparent spread in the distribution for large primary energies 
above E/Eq > 1 is due to inclusion of large zenith angles as evident from 
Figure 7. The effect of restricting the zenith angles to < 70“ is shown by 
the dotted line for Eq = 50 TeV. This effect is also noticed in Figure 8, in 
which the inclusion of angles beyond 60“ spread the distribution beyond the 
peak energy. In an actual instrument, because of the detector depth, there 
would be zenith angle restriction beyond which the instrument cannot view. It 
is clear from this figure that for energy estimation, an omnidirectional 
detector introduces no more than 15 percent dispersion towards under 
estimating the energy and no more than 5 percent towards over estimating the 
electron energy. 



The major source of uncertainty in determining the energy arises from the 
fluctuation in the synchrotron emission process* This fluctuation depends 
upon the value of n^ chosen; the larger the value of n^, the smaller is the 
fluctuation. The mean number of photons detected would be closer to the value 

of n^ for small energies and increases as the electron energy increases. The 
mean number also depends on the area of the detector. As an example, the mean 

number of photons at 10 TeV is about 6.8, for n^ “ 5 and detector radius of 
100 cms; this value increases to about 8 as the radius of the detector is 
increased to 200 cms. At 50 TeV, the corresponding values are about 9.3 and 
9.8, respectively. In order to determine the effect of fluctuation one needs 
to incorporate the detector efficiency and its energy response, and therefore 
we do not deal with this aspect here. 

In the case of a balloon-borne experiment, one may expect that low energy 
bremsstrahlung radiation by an electron in the residual atmosphere could mix 
with the synchrotron photons. However, it can be shown that the Landau- 
Pomeranchuk-Migdal effect (Stanev et al., 1982) sets in around a few hundred 
MeV, suppressing the emission of low energy photons. On the other hand, the 
effect of scattering of the electrons could increase the effective response 
area by a small amount. Although we have taken into account the attenuation 
of photons by Compton process, it may be mentioned that a small fraction of 
these scattered photons might arrive at the detector affecting the 
col linearity, when the minimum n^ criterion is applied. 

6. DETECTION OF ELECTRONS 

We have shown that the response area of the detector depends upon the 
arrival direction, n^, Eq and the active area of the detector. For an 
instrument with a finite acceptance angle, the collecting power of the 
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Instrument would depend upon ny, Eq and A. The collecting power of the 
instrument is obtained by integrating CP(Z,(j)) over Z and We have plotted 
in Figure 13, the collecting power of a detector as a function of the active 
area of the detector for various values of Oy and Eq, For this purpose we 
have considered an omnidirectional circular detector of radius R, which is 
capable of detecting photons of Ey > 20 keV, flown over Palestine, Texas at 4 
g cm"^ of atmospheric depth. Curves A, B, and C correspond to collecting 
powers calculated using Oy values of 3, 5, and 7 photons, respectively; the 
solid curves are for Eg = 10 TeV and the dashed curves are for Eg = 100 TeV. 
One notices that curve C is very much steeper than curve A indicating that the 
minimum area required to detect electrons of energy 10 TeV using ny « 7 
photons is about 1.8 m^; similarly for rty « 5 photons, it is about 0.75 m^. 

The solid curve A and all the dashed curves are parallel to each other, and 
they can be represented by a power law of the type proportional to A°‘, where a 
is about ’.43. This is an interesting feature of this technique of detecting 
electron, in which the collecting power of the instrument has a higher order 
dependence on the detector area. 

The collecting power of this detector is plotted in Figure 14 as a 
function of the energy of electrons, for various values of Oy and A. In this 
figure, curves A, B, and C correspond to different detector sizes with radii 
200 cms, 100 cms, and 50 cms, respectively. The solid and dashed curves 
represent two different selection criteria with ny = 3 and 5 photons, 
respectively. At large values of Eg, these curves become straight lines 
indicating that the collecting power is proportional to Eg. As the energy 
decreases, one can recognize the threshold effect of the detector. The 
collecting power decreases by an order of magnitude as the energy decreases by 
a factor of 2 soon after it deviates from the linear dependence in energy. 
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This behavior provides the threshold detection for electrons and, thus, the 
Integral flux can be determined accurately even by a detector with limited 
energy resolution. It Is also evident from this figure and from Figure 13, 
that the collecting power Is very sensitive to both ny and A. One Infers from 
this study that, as long as A Is not very small and ny not very large, the 
collecting power of a detector Is asymptotically proportional to A^E^nyx 
where a = 1,43 and 0 » - 1.8. 

An instrument designed to detect electrons using this technique should be 
very efficient in detecting photons above 20 keV, because the number of 
photons emitted by the electron is rather small. This can be achieved by 
assembling detector elements which respond to different energy Intervals 
efficiently. Proportional counters can be used for X-rays of energy less than 
60 keV. These can be followed by Nal or Csl crystals or scintillators to 
detect photons up to a few MeV, and they are further followed by lead-glass 
Cherenkov detectors to absorb high energy gamma-rays, Good spatial resolution 
can be achieved to any desired accuracy in the proportional chambers without 
difficulty. In the case of other detectors, one can make use of hodoscope 
techniques (for example, Streitmatter et al., 1981). One needs to make 
modules of areas less than 50 cms X 50 cms, and they are to be put together to 
achieve the required total area of the detector. Simultaneous coincidence in 
more than two such modules would eliminate the vast background of charged 
particles and reduce the chance coincidence event rates to a minimum level. 
Apart from a few events due to such chance coincidence, large angle air 
showers, gamma-ray bursts, and high energy Interactions above the detector can 
trigger the system. The electron event can be unambiguously identified from 
Its unique CO-1 Inearl ty of the observed photons. Time sequence of these 
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photons can further eliminate any chance coincidence. Such a detector can 
also be made use of for the study of gamma-ray bursts of short durations. 

For an ideal detector of the kind described above, we have calculated the 
expected rate of electron events by making use of the curves in Figure 7 and a 
differential energy spectrum of electrons. The recent data on the observed 
electron spectrum, as summarized by Stephens (1981), can be represented by a 
power law of the type jfE^) « 500 electrons/(m^ sr s GeV). The event 
rate thus calculated is shown in Figure 15 as a function of the energy of 
electron for various values of n^ and A. Curves A, B, and C in this figure 
correspond to circular detectors with radii 200 cms, 100 cms, and 50 cms, 
respectively; the solid and dashed curves relate to values of n^ = 3 photons 
and 5 photons, respectively. The collection time in a typical balloon flight 
is about 20 hours. Using a detector of area ICm^, a total of about 22 
electrons can be observed at energies above 5 TeV using n^ criterion of 3 
photons; this number reduces to about 8.5 for n^ = 5 photons. From a 14-day 
space shuttle flight one can collect about 180 events above 5 TeV and about 5 
events above lO^'^ eV. As there is no atmospher t absorption of X-rays, the 
detection threshold is decreased and the number of events would be larger. We 
have compared these event rates with that expected from a proposed shuttle 
experiment (Ormes, 1982) as shown by the curve denoted as TRIG. The study of 
electrons by the TRIG experiment is based on the principle of transition 
radiation to identify the electron and ionization calorimeter to measure its 
energy. The expected number of events from this experiment is only 0.5 
electron above 5 TeV. 
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SU MMARY AND DISCUSSIONS 

We have shown the viability of detecting electrons above a few TeV using 
a technique In which the electron can be unambiguously Identified by the 
photons It emits during Its traversal In the geomagnetic field. We summarize 
here some of the Important deductions made while examining this technique of 
detecting ultra high energy cosmic ray electrons. 

a. Synchrotron photons are emitted In a narrow cone along the 
direction of electrons so that they would be Incident almost along a straight 
line on the detector simultaneously. Thus, an electron event can be 
Identified unambiguously by the col linearity and the simultaneity of the 
photons detected. In addition, If the electron passes through the detector, 
all these photons would lie only on the convex side of the curvature of the 
electron, which makes the Identification very unique. 

b. Most of the photons have energies 1n the X- and y-ray regime. 
Therefore, one requires large area conventional X- and y-ray detectors with 
high efficiency for photon detection and spatial resolution good to about a 
few cms. The same Instrument should also be made use of for the study of 
y-ray bursts without much modification. 

c. Since the mean energy of the photons detected Is found to be 
approximately proportional to the square of the electron energy, one can 
determine the energy of the electron very accurately. However, the 
uncertainty resulting from the fluctuations on the emission process depends on 
the number of photons required to Identify the electron. It is also shown 
that. In order to determine the energy spectrum of electrons accurately over m 
decades in energy, one needs a detector capable of measuring photons over an 
energy span of (m^ + 2) decadiis. 
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d. We have demonstrated that the response of a detector is very 
sensitive to (i) the arrival direction and the energy of the electron, (ii) 
the number of photons required for the identification of an electron event, 
and (iii) the area and shape of the detector. Further, it also depends upon 
the geographic location and the altitude at which the experiment would be 
performed. Therefore, one needs to carefully evaluate the instrumental 
response, while designing a large area experiment, after knowing the detector 
characteristics and the location of the experiment. 

e. Though, in principle, one needs to know the arrival direction of 
an electron in order to determine its energy, we have shown that an 
omnidirectional detector can be effectively used for the measurement of the 
electron spectrum. The uncertainty introduced in the energy response, by an 
omnidirectional detector is asymmetric about the most probable value; the 
dispersion is about 15 percent towards the Tower energies and is about 5 
percent towards higher energies. The energy response function of such a 
detector is found to be Independent of (i) the energy of the electron, (ii) 
the detector area, and (iii) the number of photons required to identify an 
electron event. 

f. An instrument with low energy threshold for detecting photons can 
be made use of to determine the integral flux of electrons accurately, even 
with limited energy resolution. 

g. An important feature of this technique is that the collecting 
power of a detector increases with the energy of the electron. This is 
advantageous for the measurement of electrons, which have a steep energy 
spectrum. Thus, with this technique, one can cover more energy than with the 
existing techniques. 
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h* We have al-o pointed out the special advantage of large area 
detectors while using this technique} the collecting power of a detector is 
proportional to 

From the above discussions, it is quite clear that the detection of 
electrons using this technique would open up in a new dimension in the study 
of cosmic ray electrons at energies, where the existing techniques start to 
fall • 
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Note added in proof: 

"It was brought to our notice that the possibility of 
detecting high energy electrons using the synchrotron 
radiation they emit in the geomagnetic field has been 
suggested in the past by Prllutskii (1972) and licBreen 
(1977). 


Me Breen, B., Possible new methods of detecting cosmic ray 
e 1 ectrons by their synchrotron r adiation in the gepmagaetlc 
field and interstellar field, Proc. 12th ESLAB Symposium, 
ESA-SP-i24 i ‘“319 , '1977 . 

Prllutskii, 0. F., The possib ility of registering primary 
cosmic ray electrons b y means of synchrotron rad iatio n in 
the geomagnetic Piel(f, Soviet Phys. JfiTP Lett., 16, 320, 
1972. 
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CAPTIONS 


1. Energy spectra of synchrotron photons in a magnetic field = 

0.4G are shown for various electron energies. The solid curves 
are the differential spectra and the dashed ones are integral 
spectra. 

2 The normalized angular distributions of photons are shown as a 

2 

function of 8/>j;, where ij/ - me /E^, for various values of E^/E^; 
these curves are independent of the energy of electron. 

3. Angular distributions of synchrotron photons are shown. Curves A 
and A' are the differential angular distribution for values 
0.24G and 0.48G, respectively; curve C is the integral 
distribution using curve A. Curve B is the differential angular 
distribution of the energy flow for 0.24G. 

4. The upper part of the figure is an exaggerated view of the 
electron trajectory. B^Aj, B 2 A 2 ... are the tangents at points 

B 2 *.. on the trajectory of the electron. The wavy lines 
indicate the photons, which are intercepted by the detector plane 
within an angle <ti. The lower part of the figure is a schematic 
diagram describing the manner in which the response area of the 
detector is calculated (see text). 

5. The response pattern of a circular detector of radius 150 ems over 
Palestine, Texas is shown for an electron arriving at a zenith 
angle 45® at two different energies, 10 and 20 TeV. The solid 
thin curves describe the circumference of the response area for n^ 
= 5 photons and for different azimuthal directions; the dashed 
curve is for ny *= 7 photons. The thick solid curves are 
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Figure 6. 


Figure 7. 


Figure 8. 


Figure 9. 


Figure 10. 


Figure 11. 
Figure 12. 


indicative of the detector area; the curves for 20 TeV are scaled 
down by a factor of two for convenience. 

The dependence of the response area (upper part), (middle part) 
and the number of photons > 20 keV emitted per 10 km (lower part) 
are shown as a function of azimuthal angle for different zenith 
angles. The solid curves are for n^ = 5 photons and the dashed 
curve for n^ = 7 photons. 

The calculated response area Is plotted as a function of zenith 
angle for various azimuth angles as indicated by curves A, B, C, 

D, E, and F. The upper, middle parts and lower parts of the 
figure correspond to electron energies of 10 TeV, 5 TeV, and 30 
TeV, respectively. 

Histograms of the mean energy of the photons as observed by a 
detector are shown for different zenith angle restrictions for a 
circular detector of radius 150 cms with Eq = 10 TeV and = 5 
photons. The broken histograms are alternatively used for the 
sake of clarity in this figure and in the next two figures. 
Histograms of mean energy of photons are shown for different n^ 
values. These histograms include arrival directions up to 60® 
zenith; the ether parameters being the same as in Figure 8. 
Histograms of mean photon energy are shown for different circular 
detectors of radii 100 cms, 150 cms, and 200 cms as indicated; the 
other parameters being the same as in Figure 8. 

The most probable energy obtained from <E^> histograms is plotted 
as a function of the electron energy. 

Histograms of the ratio of the estimated energy of electrons using 
equation (11) to the true energy are shown for an omnidirectional 
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circular detector of radius 150 cms, with == 5 photons, for 
different values of Eq as Indicated, The dashed histogram 1s for 
the arrival directions less than 70® zenith. 

Figure 13, Collecting power of an omnidirectional detector Is plotted as a 
function of the detector area for different values of ny The 
solid curves are for Eq = 10 TeV and the dashed curves are for Eq 
= 100 TeV, 

Figure 14, Collecting power of an omnidirectional detector is plotted against 
the energy of electron for circular detectors with different 
radii. The solid curves are for n^ = 3 photons and the dashed 
curves are for n^ = 5 photons. 

Figure 15, The electron event rate expected in a detector 1s plotted against 
the electron energy, the symbols used here are the same as In 
Figure 8, The curve denoted as TRIC is for a proposed experiment 
in a shuttle using a transition radiation detector and an 
ionization calorimeter. 
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